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ABSTRACT 

We present a detailed study of the structural properties of four gravitationally lensed 
disk galaxies at z = 1. Modelling the rotation curves on sub-kpc scales we derive the 
values for the disk mass, the reference dark matter density and core radius, and the 
angular momentum per unit mass. The derived models suggest that the rotation curve 
profile and amplitude are best fit with a dark matter component similar to those of 
local spiral galaxies. The stellar component also has a similar length scale, but with 
substantially smaller masses than similarly luminous disk galaxies in the local universe. 
Comparing the average dark matter density inside the optical radius we find that the 
disk galaxies at z = 1 have larger densities (by up to a factor of ~ 7) than similar 
disk galaxies in the local Universe. Furthermore, the angular momentum per unit 
mass versus reference velocity is well matched to the local relation, suggesting that 
the angular momentum of the disk remains constant between high rcdshifts and the 
present day. Though statistically limited, these observations point towards a spirals' 
formation scenario in which stellar disks are slowly grown by the accretion of angular 
momentum conserving material. 
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Integral Field Spectroscopy, Gravitational Arcs: Individual 



1 INTRODUCTION 

It has long been known that the kinematics of spiral galaxies 
do not show Keplerian fall-off in their rotation curves, but 
rather imply the presence of an invisible mas s component 
in addition to the stellar and gaseous d isks l|Rubin et al.l 
ll98d : lBosmJll98ll : lPersic fc Saluccill 19881 ). In the local Uni- 
verse, observing the distribution of star light and mapping 
the gaseous component through Hi it has been possible to 
build up a picture of the how the baryonic component of disk 
galaxies is distributed and h ow this relates to th e underlying 
dark matter component (e.g. lPersic et al.lll996l h Tracing the 
evolution of galaxy mass from high-redshift up to the present 
day is only truly reliable if we can observe the same compo- 
nents at early ti mes. A pioneering study was performed by 
IVogt et al.lll996l where a handful of rotation curves (RC's) 
of objects at z ~ 0.5 was obtained using traditional longslit 
spectroscopy on the 10 meter Keck Telescope. However, at 
higher redshifts (e.g. z^l) galaxies are much fainter and have 
smaller angular disk scale lengths than galaxies observed at 



low redshift, therefore obtaining the spatial information re- 
quired for detailed studies is beyond the limits of current 
technology. Indeed, mapping the internal properties and dy- 
namics of both the stellar and gaseous components of galax- 
ies at high redshift is one of the main science drivers for the 
next generation of ground and spaced based telescopes at 
many wavelengths (e.g. ELT, NGST, ALMA). 

One way to overcome this problem is to use the nat- 
ural amplification caused by gravitational lensing to boost 
the size and flux of distant galaxies which serendipitously 
lie behind massive galaxy clusters. This technique is ex- 
tremely useful since we are able to target galaxies which 
would otherwise be too small and faint to ensure a suffi- 
ciently high signal-to-noise spectroscopy in conventional ob- 
servations. As such, gravitational lensing has been exten- 
sively used to make detailed studies of distant galaxies: for 
a galaxy at 0=1 with an amplification factor ten, an angular 
scale of 0.6" can correspond to ^0.5 kpc, sufficient to map 
also a small spiral. 
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The benefits of gravitationai iensing are complemented 
by Integral Field Spectroscopy (which produces a contiguous 
two dimensional velocity map at each point in the target 
galaxy). This allows a clean decoupling of the spatial and 
spectral information, thus eliminating the problems arising 
from their mixing in traditional long-slit observations. It is 
therefore much easier to identify and study galaxies with 
regular (bi-symmetric) velocity fields. 

In this paper, we present a detailed study of four 
rotation curves extracted from disk galaxies which have 
been observed through the cores of massive galaxy clus- 
ter ^Jense^JTliese targets are t aken from the recent work 
of ISwinbank et all l|2003l , 120061 ). They were observed with 
the Gemini-North Multi-Object Spectrograph Integral Field 
Unit (CMOS IFUfl 

We concentrate on the galaxy dynamics 
as traced by the [On]AA3726.1,3728.8A emission line dou- 
blet. The IFU data provide a map of the galaxy's velocity 
field in sky co-ordinates. To interpret this field the magni- 
fication and distortion caused by the gravitational Iensing 
effec t is removed using d etailed mode l s of th e cluster lenses 
(see ISmith et all l|2005l) . iKneib et all (|l996h . ISmithl (|2002| ) 
for details). The primary constraints on the lens models are 
the positions and redshifts of spectroscopically confirmed 
gravitational arcs in each cluster. The source plane velocity 
fields of four systems which display regular (bi-symmetric) 
rotational velocity fields i.e. resembling rotating disks, were 
reduced to one-dimensional rotation curves from which the 
asymptotic terminal velocity was extracted and compare d 
with the galaxy luminosity (Swinbank et al. 2003, 2006). 
The key advantage of using gravitational Iensing to boost 
the images of distant galaxies is that we are less biases to- 
wards the most luminous galaxies. Whilst observational in- 
formation on the distribution of the Hi disks in galaxies at 
these early times would be wellcome, such observations will 
have to wait for future instrumentation (e.g. ALMA). 

In this paper we use nebular emission lines to probe the 
kinematic of the galaxies. We extract one-dimensional rota- 
tion curves from the velocity fields to infer the distribution 
of stellar and dark matter components. Finally, we compare 
our results with similarly luminous disk galaxies in the local 
Universe. Through-out this paper we use a cosmology with 
H = 72kms~\ fi = 0.3 and A = 0.7, to = 13.7 Gyr. 



2 DATA 

2.1 Sample Selection 

O ur sample comes f rom o bservations of six gravitational arcs 
in lSwinbank et all |2006). In order to avoid possible biases, 
the targets were selected only to be representative of galax- 
ies in the distant Universe and no attempt was made to 
select galaxies with relaxed late-type morphology. We did, 



1 Programme ID: GN-2003A-Q-3. The CMOS observations are 
based on observations obtained at the Gemini Observatory, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under a cooperative agreement with the NSF on 
behalf of the Gemini partnership: the National Science Founda- 
tion (United States), the Particle Physics and Astronomy Re- 
search Council (United Kingdom), the National Research Council 
(Canada), CONIC YT (Chile), the Australian Research Council 
(Australia), CNPq (Brazil) and CONICET (Argentina). 
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Table 1. Derived structural parameters from the RC mass mod- 
elling. Error bars for M-q are shown in [2] while the uncertainties 
in ro, v oh , R op t = 3.2.Ro and po amount to 30%, 15%, 15%, 40% 
respectively. 



however, require that arcs were resolved in both spatial di- 
mensions so that a two dimensional velocity field could be 
extracted from the IFU data. This restricted our selection 
to galaxies with moderate magnification. From the sample 
of six galaxies, four galaxies appear to have (relaxed) bi- 
symmetric velocity fields with late-type morphologies and 
colours. The rotation curves from these four galaxies ap- 
pear regular and we therefore restrict our analysis to these 
arcs. We stress that the morphology, colours and velocity 
fields of the four galaxies in this sample all strongly sug- 
gest these galax i es are consistent with late type spirals (see 
ISwinbank etail i|2006h ). 



2.2 Photometry 

From our optical/near- infrared imaging, we constrain the 
spectral energy distribution (SED) of each galaxy. Since the 
arcs usually lie with a few arcs-seconds of nearby bright 
cluster galaxies, we calculate the magnitude of the arcs in 
various pass-bands by masking the arc and interpolating the 
light from the nearby cluster members. The background light 
is then removed and surface photometry in different bands 
is obtained. 

Using the cluster mass models the arcs are recon- 
structed to the source-plane and the geometry and disk-scale 
parameters of the disks are measured. This is achieved by 
fitting ellipses to an isophote of the galaxy image using a 
modified version of the idl GAUSS2dfit routine which fits 
an exponential profile to the two-dimensional light distribu- 
tion. From this, the ellipticity, the inclination and luminosi- 
ties and the disk scale leng ths are obtained (see Table 2 and 
3 in ISwinbank et af] (|2006t )). These latter are also reported 
below in Table 1. 



2.3 One-Dimensional Rotation Curves 

In Fig. [1] we show the one dimensional rotation curves of the 
galaxies in our sample. These are extracted by sampling the 
velocity field along the major axis cross section. The zero- 
point in the velocity is defined using the center of the galaxy 
in the reconstructed source plane image. The error bars for 
the velocities are derived from the formal 3<r uncertainty in 
the velocity arising from Gaussian profile fits to the [On] 
emission in each averaged pixel of the datacube. For the 
mass modelling analysis we folded the rotation curves on 
the kinematical center to ensure that any small-scale kine- 
matical sub-structure is removed. 
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Figure 1. Filled circles represent the IFU rotation curves having enough data to allow the coaddition of the kinematics of the receding 
and the approaching arm, red and blue open circles represent the rotation curves along each arm in the other cases. The dark-matter 
and stellar components are shown with long dash and short dash lines, whilst the model circular velocity is shown with a solid line. 



3 MODELLING AND RESULTS 

The rotation curves of local spiral galaxies imply the pres- 
ence of an invisible mass component, in addition to the stel- 
lar and gaseous disks. The paradigm is that the circular 
velocity field can be characterized by: 

V 2 = vl + yf + y| z (1) 

where the subscripts denote the stellar disk, dark halo 
and gaseous disk respectively. From the photometry we 
mod el the stellar c omponent with a Freeman surface den- 
sity (|Freemanlll970T ): 

where Rd is the disk lenght-scale, while R op t = 3.2Rd can 
be taken as the "size" af the stellar disk, whose contribution 
to the circular velocity is: 

Vg(x) = i^(3.2x) 2 (/ ()J rY - hKi) (3) 

where x = R/R op t = R/(3.2Rd) and I n and K„ are the 
modified Bessel functions computed at 1.6 x. In small spirals 
the Hi disk is an important baryonic component only for R > 
3Rd, i.e. outside the region considered here. It is plausible 
that in similar objects also at high redshifts, inside 3Rd, 
the Hi gas contributes to the gravitating baryonic mass by 
a very small amount, residing at larger radii and it slowly 
infalls forming the stellar disk. 

For the dark matter component we take a spherical 
halo for which V„{R) = G M H (< R)/R. Following the 
observational scenario constructed in the local Universe 



|Persic fc Saluccilll988l; ISalucci fc Burkerdlioool ) we assume 
that it has the iBurkertT i 1995T ) density profile (see also 
Salucci et al 2007): 

p(R) = 



Po ro 



where ro is the core radius and po the effective core density. 
It follows that: 



M H (R) = k 



+ 2 ln 



(R- 



- ro) (R 2 + rg) 



(4) 



(5) 

with k = GAp rl and of course Vg(R) = GM H (R)/R We 
note that the adopted velocity profile is a quite general: it 
allows a distribution with a core of size ro, converges to the 
NFW profile at large distances and, for suitable values of ro 
can also mimic the NFW or a isothermal profile, over the 
limited region of galaxy which is mapped by the rotation 
curves. 

The mass model has three free parameters: the disk 
mass Md, the core radius ro, and the central core density po- 
The observations extend out approximately (2 — 3)Rd, have 
10-60 independent measurements with an observational er- 
ror of 3% — 10% in their amplitude, of 0.05-0.2 in their slopes 
dlogV / dlogR. The error in the estimate of the disk inclina- 
tion angles is negligible with respect to the above. These er- 
rors are (understandably) higher than those associated with 
the best quality local RC's and make it difficult to constrain 
the halo density profile p(R). However, they are sufficiently 
small to yield a reliable value for the halo mass, the aver- 
age density inside a reference radius, (which we chose to be 
R 2 = 2R D so that < p >= M H (7? 2 )/(4/37r7?i) and the disk 
mass and a reasonable estimate of the " core radius" . 

By reproducing the observed rotation curves with the 
models given by equations 1-5 we derive the best fit param- 
eters for each galaxy and overlay the resulting mass model 
onto each rotation curve in Fig. [T] In Table 1 we report the 
main structural parameters: the disk mass, the halo core ra- 
dius, po, the optical radius R op t and voh = VniRopt), i.e. the 
halo contribution to circular velocity at the optical radius. 

In all our rotation curves the amplitude and the pro- 
file of the stellar disk contribution can not reproduce the 
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Figure 2. Left: disk mass (in solar units) versus reference velocity (in km/s) compared with the 2=0 relation (solid line) and with this 
scaled down by 0.6 dex (dashed line). Middle: average dark matter density (in g cm -3 ) versus reference velocity compared with the 2=0 
relation (solid line). Right: disk angular momentum per unit mass j (in km/s kpc units) versus disk mass compared with the 2=0 relation 
(solid line). The uncertainties on the local relations are: 0.15 dex in j, 0.2 dex in log Mjj and 0.3 dex in log < p >. 



observed RC rise between 1.5 Rd and the last measured 
radius. This strongly suggests evidence for the presence, at 
z ~ 1, of a dark matter component of mass comparable 
to that found for local disk galaxies with similar V op t (see 
Figures 2, 8 and 9 of Persic, Salucci & Stel, 1996). We de- 
rive disk masses ranging between 1 x 1O 9 M to 2 x 10 10 M Q 
for galaxies with reference velocity V op t between 100 km s -1 
and 200kms _1 . These disk masses are smaller by a factor 
2-4 than those of the local spirals with the same reference 
velocity which are shown to follow Inner Baryon Dominance 
(see Salucci & Persic, 1999). The high redshift stellar disks 
are sub-maximal disks. Forcing a maximal disk even in the 
"weak" implementation of Persic and Salucci (1990) leads 
to unacceptable fits of our high-z RC's. 

The best fit values for ro is of the order of ~ 1.5-Rd, 
which is larger than usually compatible with a NFW pro- 
file, although the error-bars on our data preclude any strong 
statement. 

Since the gravitational lens model affects the source- 
plane reconstruction of the galaxies, we must test for the 
effect that this has on the resulting rotation curves. We 
reconstruct the galaxies using the family of lens models 
that inhabit the A\ 2 contour corresponding to la confi- 
dence interval relevant to each cluster. For example, the 
model of RGB 1745 has five free parameters, and the lens 
modeling uncertainties are therefore derived by consider- 
ing models within the Ax 2 =5.89 contour. From each re- 
construction we extract the one-dimensional rotation curve 
and apply the analysis outlined above and find the maxi- 
mum variations are: A log Mb = 0.03, A log < p >= 0.02, 
Alog(V" op i) g 5kms . Thus the uncertainties in the gravi- 
tational lens modelling is negligable compared to the uncer- 
tainties in the RC mass modelling. 

A cosmological significance of our result is evident in 
Fig[2]where we compare the disk mass, the mean dark matter 
density within the optical radius and the angular momen- 
tum per unit mass, all as a function of reference velocity. 
These are compared to si milar properties of the local ob- 
jects (IShankar et al.ll2006l ). Figure [2] strikingly shows that 
high redshift spirals, modulo an offset of 0.6lg 15 dex, are 
on the same logMo versus logV op t relationship found for 
local spirals arising from the systematic structural proper- 



ties of their mass d istribution l|Tonini et al.|[2006T ). see also 
l|Salucci et al.lll993l ) . 

In Fig 2, from the values of < p >, a quantity that 
differently from po is weakly affected by the RC \-a fitting 
uncertainties, it is apparent that the DM halos of z — 1 disk 
galaxies are denser by O.Tto'J dex than those around sim- 
ilarly luminous z=0 spirals. The evidence that spiral disks 
at z=0 and z=l have the same structural relationships is 
further supported by observations of the evolution of the 
Tully-Fisher relation (which correlates the disk mass with 
Vopt)- In our and in other independent samples (Swinbank 
et al 2006, , Vogt et al 1996, Bamford et al 2005) the galaxies 
at z = 1 show a TF relation with a slope similar to that of 
the local TF, but with an offset compatible with that found 
in the present work from the disk mass vs rotation velocity. 
This suggests that from z=l to z=0, the stellar disk masses 
Mo of a spiral has grown by a factor ~ 4l 2 > that leads to 
just a modest increase in the DM dominated quantity V op t ■ 

Further evidence that z = 1 disk galaxies are related to 
present day spirals is provided by the relationship between 
angular momentum per unit mass (j) versus the reference 
velocity, as shown i n Fig. [2] This wel l theoretically moti- 
vated relation (e.g. (|Tonini et al.| [2006;l) can be considered 
as the imprint of the process of the formation of disks in- 
side dark matter halos related to the cosmological proper- 
ties of halo spin parameters. As Figure [5] shows, there ap- 
pears to be no evolution in this crucial relationship between 
the the comological time at which we observe these spirals, 
2 = 1,4 = 6 Gyr and the present time, z = 0, t = 13.7 Gyr. 
This agreement is remarkable: it establishes a link between 
local and high redshift disks, supporting the idea that the 
angular momentum remains constant during the evolution 
of a disk system from high redshift to the present day. 

The kinematical estimate of the disk mass allows us to 
derive the mass-to- light ratios for our disk systems as a func- 
tion of luminosity and colour. In Fig. 3 we compare the disk 
mass and colour as a function of ma ss-to-light ratio com - 
pared to the relation in local spirals (|Shankar et al.l 12006). 
In order to compare directly with local relations, we consider 
a simple passive evolution model for the luminosity evolu- 
tion. For a single stellar population the zero-point of the 
local relation is decreaed by a factor log(13.7/6) which ac- 
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Figure 3. Left: Mass-to-light ratio versus disk mass for the galaxies in our sample compared with the local disk galaxies from Shankar et 
al. (1996) (solid line). The average offset corresponds to an excess in luminosity by a factor four. Right: Mass-to-light ratio versus colors 
relations for our sample. The square box denotes the predicted colours from a Bruzual & Chariot single stellar population (see text for 
details). 



counts for the passive evolution of a single stellar population 
from z = 1 to z — 0. As Fig. 3 shows the mass-to-light ratios 
as a function of galaxy (B — V) colour are in broad agree- 
ment with p redictions of a single stel lar population which is 
~ lGyr old IjBruzual fc Charlotll2003l ). although clearly pho- 
tometry at other wavelengths (such as rest-frame K-band) 
would allow a more detailed decomposition of the stellar 
popultations in these galxies. 



3(M*(t )/(W w MQ)M Q /yr which points towards a declin- 
ing SFR history. 

The marked increase of the luminosity per unit stellar 
mass in objects at high redshifts with respect to their local 
counterparts has the simplest explanation in a passive evolu- 
tion of the starforming disks. Obviously, this simple picture 
requires us to assume that the high redshift systems are the 
direct counter-parts of similar rotation speed spirals at low 
redshift. 



3.1 Cosmological Implications 

Figure 2 shows that at a radius corresponding to V op t the 
high redshift galaxies are significantly denser than compara- 
bly luminous local disk-galaxies: the average offset is about 
0.6 dex in log(< p) >. Although we can not exclude that 
dynamical processes occur between z = 1 and z — to re- 
duce the dark-matter density in the luminous regions, this 
offset is naturally explained if the halos embedding these 
disk galaxies formed at earlier times than the halos around 
similarly massive z = spirals. In this framework we esti- 
mate the ratio between the virialization redshift of the local 
galaxies and and that of the galaxies in our sample. Since 
p v oc A(z v )(l + z v ) 3 where p v and z v are average density 
and redshift at virialization and A„ is known, for zo = 1, 
z v = 1.7, which corresponds to t v — 6 Gyr. 

Assuming that our sample is a fair representation of 
disk galaxies at z ~ 1 and that these are approximately co- 
eval, from the comparison of their structural properties with 
those of z — spirals, the following simple picture emerges: 
a present day spiral, with a given circular velocity, half-light 
radius and the angular momentum per unit mass, at redshift 
1 had similar values for these quantities, but a smaller stellar 
mass: < M+(t i, a ) /M*(to) >— 0.3. This induces a scale for 
the average SFR in the past 8 Gyr: ~ 0.75AL(t )/ (t obs - 
t ) ~ l(Af*(t o )/(lO lo M )M Q /yr. With these disks hav- 
ing an average age of 1 Gyr at z — 1 we can also derive 
an "early times" average SFR ~ 0.25M* (t )/(lGyr) ~ 



4 DISCUSSION AND CONCLUSIONS 

In this study, we have investigated the detailed properties 
of four disk galaxies at z = 1. These galaxies were observed 
at high spatial resolution thanks to the boost in angular 
size provided by gravitational lensing by foreground massive 
galaxy clusters and allow a much more detailed comparison 
with local populations than usually possible for galaxies at 
these early times. M odelling the one-di mensional rotation 
curves with those of iPersic et al. I (|l996h we derive best fit 
parameters for the total dynamical mass, the core radius, 
the effective core density and the angular momentum per 
unit mass. 

The best fit model rotation curves to the data show that 
the amplitude and profile of the stellar disk componentcan 
not unambiguously reproduce the rise in the rotation curve 
without a dark matter component. Comparing the average 
dark matter density inside the optical radius we find that 
the disk galaxies at z = 1 have larger densities (by up to 
a factor of ~ 7) than similar disk galaxies in the local Uni- 
verse. In comparison, we find that the angular momentum 
per unit mass versus reference velocity is well matched to 
the local relation suggesting that the angular momentum 
of the disk remains constant between high redshift and the 
present day. Though statistically limited, these observations 
point towards a spirals' formation scenario in which stellar 
disk are slowly grown by the accretion of angular momen- 
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turn conserving material. Our result, also consistent with 
the theoretical evolution of the angular momentum of disks 
from semi-analytic models from z=l to z=0 which show the 
modest offset of Aj^0.2 kpc. kms -1 for objects with circu - 
lar velocities betwe en 50 and 300 kms -1 |Cole et alj l|2000(k 
iBower et al.l ([2006) ) . These reults provide an evolutionary 
link between the disk systems we observe at redshift z ~1 
and the present day population of spirals. 

While these results are based on data of only four ob- 
jects, they nevertheless show the power which gravitational 
lensing can have on studying the internal properties of high 
redshift galaxies. In particular, these observations should be 
viewed as pathfinder science which will soon be routine with 
Adaptive Optics Integral Field Unit observations on eight 
and ten meter telescopes (e.g. iGenzel et alj [20061 ). More- 
over, by combining these results with upcoming telescopes 
and instruments (e.g. ALMA) which will be sensitive enough 
to map the Hi content of galaxies to z = 1, the exact relation 
between gas, stars and dark matter can be probed in much 
more detail. 



Shankar, P., Lapi, A., Salucci, P., De Zotti, G., & Danese, 

L. 2006, ApJ, 643, 14 
Smith, G. P. 2002, Ph.D. Thesis, University of Durham, 

UK 

Smith, G. P., Kneib, J.-P., Smail, L, Mazzotta, P., Ebeling, 

H., & Czoske, O. 2005, MNRAS, 359, 417 
Swinbank, A. M., Bower, R. G., Smith, G. P., Smail, I., 

Kneib, J.-P, Ellis, R. S., Stark, D. P., & Bunker, A. J. 

2006, MNRAS, 368, 1631 
Swinbank, A. M., Smith, J., Bower, R. G., Bunker, A., 

Smail, L, Ellis, R. S., Smith, G. P., & Kneib, J. P., et al. 

2003, ApJ, 598, 162 
Tonini, C., Lapi, A., Shankar, F., & Salucci, P. 2006, ApJL, 

638, L13 

Vogt, N. P., Forbes, D. A., Phillips, A. C, Gronwall, O, 
Faber, S. M., Illingworth, G. D., & Koo, D. C. 1996, ApJL, 
465, L15 



ACKNOWLEDGMENTS 

We would like to thanks the anonymous referee for his/her 
suggestions which improved the content and clarity of this 
paper. We thank Carlton Baugh for providing the theoretical 
evolution of disk mass in galaxies from GALFORM. We also 
thank Gigi Danese for useful discussions. AMS acknowledges 
support from a PPARC Fellowship, RGB acknowledges a 
PPARC Senior Fellowship and IRS and GPS acknowledges 
support from the Royal Society. 



REFERENCES 

Bamford, S. P., Milvang- Jensen, B., Aragon-Salamanca, 

A., & Simard, L. 2005, MNRAS, 361, 109 
Bosma, A. 1981, AJ, 86, 1825 

Bower, R. G., Benson, A. J., Malbon, R., Helly, J. O, 
Frenk, C. S., Baugh, C. M., Cole, S., & Lacey, C. G. 2006, 
MNRAS, 659 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Burkert, A. 1995, ApJL, 447, L25 

Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, 

MNRAS, 319, 168 
Freeman, K. C. 1970, ApJ, 160, 811 

Genzel, R., Tacconi, L. J., Eisenhauer, F., Forster 

Schreiber, N. M., Cimatti, A., Daddi, E., Bouche, N., & 

Davies, R., et al. 2006, Nature, 442, 786 
Kneib, J.-P, Ellis, R. S., Smail, I., Couch, W. J., & 

Sharpies, R. M. 1996, ApJ, 471, 643 
Persic, M. & Salucci, P. 1988, MNRAS, 234, 131 
Persic, M., & Salucci, P. 1990, MNRAS, 247, 349 
Persic, M., Salucci, P., & Stel, F. 1996, MNRAS, 281, 27 
Rubin, V. O, Thonnard, N., & Ford, Jr., W. K. 1980, ApJ, 

238, 471 

Salucci, P., Frenk, C. S., & Persic, M. 1993, MNRAS, 262, 
392 

Salucci, P. & Burkert, A. 2000, ApJL, 537, L9 
Salucci, P., & Persic, M. 1999, A&AP, 351, 442 



